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Abstract 

We report exact numerical diagonalization results of the infinite-range Ising 
spin glass in a transverse field T at zero temperature. Eigenvalues and eigen- 
vectors are determined for various strengths of F and for system sizes N < 16. 
We obtain the moments of the distribution of the spin- glass order parameter, 
the spin-glass susceptibility and the mass gap at different values of T. The 
disorder averaging is done typically over 1000 configurations. Our finite size 
scaling analysis indicates a spin glass transition at Fc — 1.5. Our estimates 
for the exponents at the transition are in agreement with those known from 
other approaches. For the dynamic exponent, we get z = 2.1 ±0.1 which is in 
contradiction with a recent estimate (z = 4). Our cumulant analysis indicates 
the existence of a replica symmetric spin glass phase for F < Fc. 

PACS nos : 75.10.Jm, 75.10 Nr, 64.60.Cn, 64.60.Fr 







Quantum phase transitions in disordered systems have been studied intensively in recent 
years [ij. Of particular interest is the Ising spin glass model in a transverse field which 
provides a rather simple model where it can be shown that at zero temperature, the spin- 
glass ordered ground state is destabilised by quantum fiuctuations. As a result, if one varies 
the strength F of the transverse field, there can be a phase transition between the spin-glass 
ordered and disordered ground states at a critical value of F = Fc. The case of infinite 
range model is specially interesting since in the absence of the transverse field, the model 
reduces to the usual Sherrington-Kirkpatrick (SK) model 0. The classical SK-model is very 
well studied. It has a finite temperature transition to a low temperature spin glass phase 
where replica symmetry is broken Q]. In presence of the transverse field F, the spin glass 
ordering occurs at lower temperatures and above a critical value of the field, the quantum 
fiuctuations destroy the order. At non-zero temperatures, however, quantum fiuctuations 
are unimportant in determining the critical behavior. In this letter, we are interested in the 
SK-model in a transverse field at zero temperature, where the transition between the spin 
glass to disordered ground state is driven by quantum fiuctuations alone through the tuning 
of the strength of the transverse field F. 

There has been quite a number of studies on this model at recent times. The perturbation 
expansion of the free energy by Ishii and Yamamoto gives the value of the critical field 
F ~ 1.5 and the susceptibility exponent ~ 0.5. The non-perturbative analysis the mean- 
field theory of quantum rotors [0 and the numerical techniques such as quantum Monte 
Carlo method give consistent estimates of the critical field and the critical exponents 
except the value of the dynamic exponent z. The Monte Carlo study suggests the value of 
z = 4, much too higher than the value z = 2 obtained in other studies. In a recent work, 
the Schrodinger equation for the model has been solved numerically and the interaction 
energy and longitudinal susceptibility are determined. The estimates of the critical field 
and the exponents from this study match with the analytical results. But, in this method 
the exponents are obtained not directly individually, but rather through a scaling relation. 
An important observation that comes out of the mean field theory is the possibility of a 
replica symmetric spin glass phase in quantum transition below F << F,, [l^J^l- This point 
certainly deserves further attention. It may be noted that most of these results are obtained 
using the approximate classical mapping of the quantum systems. 

Here we report the results of exact diagonalization of the infinite range Ising spin glass 
model in transverse field, using the Lanczos algorithm for system sizes < 16. Config- 
uration averaging is done typically over 1000 realizations of bond values (for low A^). We 
obtain the various moments of the order parameter distribution and the mass gap. 

The cumulant analysis of Pifl) provides an estimate of the critical field F^ — 1.5 in the limit 
of large A^. This value of Fc is comparable to the values obtained from other studies. For 
spin glass susceptibility and the mass gap, our finite size scaling analysis follows the idea 
of the phenomenological renormalization group transformation []lT|. Our estimate of the 
correlation length exponent v = 0.252 ± 0.004, the susceptibility exponent 7 = 0.51 ± 0.02 
and the dynamic exponent z = 2.1 ±0.1 (obtained from the scaling of the mass gap |r2[ ) 
compares well with the values obtained from other studies. Our estimate of z does not 
match with the quantum Monte Carlo result ||^ which we think is too high and may be an 
artifact of the size effect along the Trotter direction. Our cumulant analysis shows that for 
F << Fc, P{q) attains a two peak structure (related by the spin inversion symmetry), the 
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peaks becoming narrower with A^. This is similar to what one gets in the classical SK-model 
above the Almeida-Thouless line (see 0]) where the glass phase is replica symmetric. The 
extended width of P{q), characteristic of replica broken phase, is clearly absent in our study. 
The SK model in a transverse field is described by the Hamiltonian 

N N 

H=Y. J.,sts^-vY.s^ (1) 

<ij> i 

where the Jjj's are long ranged and follow the Gaussian distribution 

Z)(j ) = (^)i/2exp(— 1^) 

and < ij > in the summation denotes that each pair of spins is taken only once. S'f and Sf 
are Pauli spin matrices. 

We have performed exact diagonalisation of the Hamiltonian (1) using Lanczos algorithm 
for = 4, 6, 8, 10, 12, 14 and 16 and repeated for various realisations of Jij values. For each 
realization of J^, all the relevant quantities like the order parameter etc. are computed and 
then averaged over all the realizations. For the diagonalisation precedure, the basis states 
1 0a > are chosen to be the eigenstates of the spin operators {S-,i = 1, N}. In the absence of 
the transverse field term, the Hamiltonian (1) reduces to the Hamiltonian of the S-K model 
and is diagonal in this representation. With the transverse field present, each eigenstate of 
(1) is obtained as a superposition of the basis states which are 2^ in number (corresponding 
to S'f = ±1). The nth eigenstate for the quantum Hamiltonian is written as 

2iv 
a=l 

and the ground state is denoted by 1-00 >• 

We calculate the moments of the distribution P{q) of the order parameter q = 
{l/N)J2i< Sf >^. The overhead bar indicates the configuration average and the < ... > 
denotes the expectation values at the ground state. P{q) is not directly obtained here, rather 
the moments rrik = q^ = q^P{q)dq are calculated using ^, 

■1 N N 

^^ = ]^E--E<^n-^^>^- (2) 

It can be easily shown that the RHS of the above equation gives the /cth moment of the 
distribution. We take for example 
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^2 = T7^EE[<^n^.^><^n5.^>] 
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Now, < SI_^Sf^ > in the ground state is given by < ipo\Sf_^S-j4'o >. The latter quantity 
is again given in terms of the basis states, i.e., < ipolSf^SfJi/jQ >= J2'^ Ic^aP < ^f^Sf^ >q,, 
where < ... >a denotes the product of the iith and 22th spin in the ath basis state. Writing 
|a° p = Ua, the RHS of (2) is therefore given by 
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■1 iV Af 2^ 2^ 

After configurational averaging, it gives m2 = J dqPj{q)q^. 

The exact diagonalisation procedure gives us the cUq's. Using the above, we have calcu- 
lated the Binder cumulant |jl3| for a value of and F as 

and the nonlinear susceptibility 

N 

X = (1/iV) Y: < S^S] . (4) 

We calculate the order parameter and the spin glass susceptibility which are given by the 
first and second moments of the distribution P{q). While calculating the moments of q, we 
take a set of half of the basis states, ignoring the other half, which consists of the trivially 
degenerate ones and obtained by simply flipping the spins in the states of this set. This is 
necessary as otherwise, even in the ferromagnetic case one will get < 5*/ >= 0, as the two 
degenerate states are equally probable. The quantities have to be suitably rescaled to get 
the proper values. 

The energy eigenvalues of the nth state is given by En =< ipnlHlipn >■ The mass gap is 

AE = El- Eo (4) 

where Eq is the ground state energy and Ei is the energy of the first excited state. We 
determine the average mass gap Ai^^7v(r) for different values of of F and N. 

We estimate first the cooperative energy per spin, which is the expectation value of the 
cooperative part of the Hamiltonian in the ground state (in this case, it is simply obtained 
by operating the first term of the Hamiltonian (1) on {ipQ >). The asymptotic value of this 
quantity, which gives the classical ground state energy per spin in the S-K model (in the 
limit r — > 0), is found to be equal to —0.74 ± 0.01 (in units of J). This agrees very well 
with the ground state energy (—0.76) of the classical SK model 0. 

The Binder cumulant gN{^) as a function of F is shown in Fig. 1 for different values of 
A^. The form of 5'Ar(F) changes with the system size and near Fc is expected to scale with 
A^ in the following way 

g^iT) ~ /((F - FJAT-) (5) 

where f{X) is a scaling function and x is related to the mean field correlation length exponent 
u. 5'7v(F) versus F for various A^ then yield a family of curves. These curves intersect at a 
common point which gives the value of Fc. Due to the corrections to finite size scaling, the 
intersections of the QnO^) curves for all A^ values may not coincide at a common point (see 
0] for discussion). We determine the point of intersection Fc(A^, A^') of the g(T) curves for 
two system sizes A^ and A^' and plot it against 1/(A^A^')^^^ in Fig. 2. The value of Te{N, N') 
converges to F^ ~ 1.5 ± 0.1 in the limit of N, N' ^ oo. This value of F^ is comparable to 
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the value obtained from the non-perturbative analysis [^], the Pade treatment and the 
solution of the Schrodinger equation 

The finite size scaling concept of Fisher and Barber extended to the long-ranged systems 
suggests that the exponent x = l/i^Umfdc), where z/^/ is the mean field correlation length 
exponent and dc is the upper critical dimensionality of the corresponding short range system 
||T5| . For any pair of systems with total spins and A^', we first find out the point of 
intersection Tc{N,N'). With Fc = Tc{N,N') we try to scale 5'Ar(g) and gN'ig) around Fc 
with a suitable value of x, so that all the data points fall on the same scaling curve f{X). 
Fig. 3 shows the value of gainst l/y/WW. The value of X is almost independent of N 
and we get x = 0.66 ± 0.01. If we take dc = Q which is the upper critical dimension of the 
classical Ising spin glass we get u = Umf = 0.252 ± 0.004 in agreement with the predicted 
value of z/ = 1/4 

The finite size scaling form for the spin glass susceptibility x can be written as 

X^(F)~iV^$((F-F,)iV-) (6) 

where the exponent y = xj, 7 being the exponent which characterizes the divergence of 
the susceptibility at Fc in the thermodynamic limit. is a scaling function. The finite 

size analysis can be done invoking the idea of the phenomenological renormalization group 
transformation (see |1T[ for a review). If is a power function of X, then the above 



scaling form is satisfied exactly for any two finite systems of spin numbers and A^' in the 
sense that Xiv(r)/A^^ = XN'iT')/N'y with the recursion relation F - F^ = (F' - Fc)(A^VA^)' 
The fixed point F* of the transformation satisfies 



XN'iT*) \N'J ' 

where F* — F^ as A^, A^' — >■ 00. For any two sizes A^ and A^', we take F* from the intersection 
of the curves QNiX) and gN'(X) and apply the above fixed point equation to get the value 
of y. The resulting estimate of y depends on A^ and A^' and approaches the exact value as 
N,N' ^ 00 (see [jT6[ for discussion). Fig. 4 shows the value of y plotted against l/Nm, 
where A^^ = \flNN'), from which we determine the asymptotic value of y = 0.33 ± 0.01. 
From y and x, we get 7 = 0.51 which agrees with the value 1/2 obtained in other studies 

m- 



The scaling form for the mass gap AE]\f{T) can be written as |T2| 

AEn{T) ~ A^^-"^((F - FJA^^) 

where \1/ is a scaling function. Employing the same scaling analysis as above, for the mass 
gap, we obtain the value of the exponent z. Fig. 5 shows the value of z plotted against 
1/Nfn- We find the asymptotic value of value 2; = 2.1 ± 0.1. This value of z agrees with its 
predicted value in |^ but disagrees with the quantum Monte Carlo simulation result z = 4. 

Proceeding in the same way, we find jS = 1.0 ± 0.1 where f5 is the order parameter 
exponent, i.e., in the thermodynamic limit g ~ (F — Fc)^- This is in agreement with the 
previous studies 

/,From Fig. 1, we find that the behavior of the Binder cumulant is very similar to what 
happens in normal spin system. If g{T) is the value of the cumulant in the thermodynamic 
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limit at any temperature T, then for a normal ferromagnetic system it is expected that 
g(T) = 1 for T << Tc and g(T) = for T >> T^. The value of g(Tc) at Tc depends 
on the dimensionality of the system. At low T-values, the value attained by the Binder 
cumulant depends on N and extrapolates to the value unity in the thermodynamic limit. 
This coincides with the fact that the order parameter distribution (P(m) of magnetisation 
m for a ferromagnet) reduces to a delta function (modulo symmetric operation) in the limit 
N ^ oo and the system goes to a definite symmetry broken state. 

This is not the case where replica symmetry is broken like in classical SK model at 
low temperatures and zero external field or spin glass models in higher dimensions. Here 
P{q) does not reduce to a delta function in the ordered phase even at ^ oo. Instead, 
it consists of a delta function plus a continuous, almost size independent part going right 
down to g = 0. As a result, in the spin-glass phase, the Binder cumulant attains a value 
which is significantly lower than unity and does not show any change with the system size 

Our g{r) curves in Fig. 1 for different A^ values show a crossing point (corresponding to 
a phase transition) at a finite value of F = Fc, below which the curves splay out for different 
A^ and then saturate to values which tend to unity with increasing A^. This suggests that for 
F << Fc, the system goes to a definite replica symmetric spin-glass ordered ground state. 

In summary, all our estimates for the critical tunnelling field (critical point) and expo- 
nents agree with the previous estimates, our estimate for z agrees with that estimated by 
Miller and Huse and also Ye et al 0], while it disagrees that obtained by Alvarez and 
Ritort [§. We believe, the high value of z obtained in the quantum Monte Carlo study (of 
Ritort et al) is due to an artifact of the size anisotropy in the Trotter direction. We also 
emphasize that the variation of gN{^) functions with A^ for small F indicates the quantum 
spin glass phase (for F below Fc) is replica symmetric in this model. 

We are grateful to M. Acharyya, S. M. Bhattacharjee, A. Dutta, H. Rieger and R. R. dos 
Santos for many discussions and suggestions. We also thank D. Dhar for a careful reading 
of the manuscriot and comments. PS acknowledges support from SFB 341. 
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FIGURES 

FIG. 1. The variation of 5'Ar(r) is shown against the transverse field T for different system 
sizes N. Note that the intersection of the curves shift towards larger T values as N is increased. 

FIG. 2. The effective critical transverse field TciN,N + 2) against 1/Nm = l/^/N{N + 2) are 
shown. Fc for the infinite system is found to be at 1.5. The best fit line is shown. 



FIG. 3. The effective values of x against l/N^ = 1/\^N{N + 2) are shown to be fairly inde- 
pendent of the system sizes. 

FIG. 4. The effective values of y arc shown against l/Nm = l/^N{N + 2). The best fit curve 
with an asymptotic value y = 0.33 is also shown. 



FIG. 5. The effective values of z against 1/Nm = 1/ \/N{N + 2) are shown. The best fit curve 
with an asymptotic value z = 2.1 is also shown. 
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